Faraday spectroscopy of atoms confined in a dark optical trap 
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We demonstrate Faraday spectroscopy with high duty cycle and samphng rate using atoms con- 
fined to a blue-detuned optical trap. Our trap consists of a crossed pair of high- charge- number 
hollow laser beams, which forms a dark, box-like potential. We have used this to measure transient 
magnetic fields in a 500/im-diameter spot over a 400 ms time window with nearly unit duty cycle 
at a 500 Hz sampling rate. We use these measurements to quantify and compensate time- varying 
magnetic fields to ~10 nT per time sample. 

PACS numbers: 33.55. +b, 07.55.Ge, 37.10.Gh 



I. INTRODUCTION 

A spin-polarized atom sample is strongly birefringent for 
near-resonant light. This magneto-optic polarization ro- 
tation can be used for sensitive alkali-vapor magnetome- 
try [1 , and has been the subject of several recent studies 
in a variety of cold atom samples [2j [3l HI |5l [6] . When 
applied to localized cold atom ensembles, the result can 
be sensitive magnetometry with linear spatial resolution 
of a few tens of microns [7 . These magnetic micro- 
scopes could be of use for imaging fields near a variety 
of surfaces, including integrated circuits [H [9] and atom 
chips designed for cold atom interferometry [10]. At a 
more fundamental level, Faraday spectroscopy has also 
been considered for searches of atomic electric dipole mo- 
ments (EDM) [11 and for nondestructive quantum state 
estimation and preparation [6l [121 [12 • Such measure- 
ments benefit from large atom numbers, long interroga- 
tion times, and "field-free" confinement, i.e. confinement 
in which the trapping potential minimally perturbs the 
measurement. 

A simple way to achieve field-free conditions for a cold 
atom sample is to release the atoms from a trap and 
probe them during freefall. A drawback of this is that 
the maximum interrogation time is limited to a few tens 
of milliseconds as the atom cloud falls away from the in- 
teraction region. Isayama et. al. [2 reported a Faraday 
signal from atoms in freefall with a 1/e decay time of 11 
ms. This limitation has been overcome by confining the 
atoms to the antinodes of a red-detuned optical lattice 
in which one of the lattice beams also serves as a probe 
beam When the atoms were held in the intensity 
nodes of a blue-detuned lattice, dark-field confinement 
was achieved, although the signal was reduced because 
the interaction with the probe was correspondingly di- 
minished. 

By confining atoms in a blue-detuned trap, however, 
it is possible to achieve the simultaneous conditions of 
long interrogation time, low-field confinement, and large 
atom- number [Mj [151 HSl [13 • Blue-detuned traps pro- 



duce lower light shifts and photon scattering rates than 
red-detuned traps, enabling deep, large volume traps 
with low power requirements. Although these traps have 
been proposed for use in magnetometry ^ [2] and EDM 
searches ^Tj, to the best of our knowledge, no experi- 
mental demonstrations have been performed. 

In this paper, we report the use of dark optical traps 
to confine atoms in a submillimeter, box-like volume for 
dynamic magnetometry using Faraday spectroscopy. The 
traps are formed from crossed, high-charge-number hol- 
low laser beams [18, 19 . By repetitively spin-polarizing 
the confined sample, we extend the measurement time 
from only a few milliseconds to ms in a single load- 

ing cycle with up to 1 kHz sampling rate. We demon- 
strate the technique by measuring and compensating am- 
bient time- varying magnetic fields, such as those arising 
from eddy currents and the AC power line. We also show 
that nonlinear spin dynamics due to the probe beam |T2] 
are preserved in these traps. The increase in duty cy- 
cle demonstrated here is promising for both magnetome- 
try and for efficient quantum state preparation based on 
these nonlinear dynamics. 




FIG. 1: (Color online) Layout of experiment. Crossed hollow 
beams confine atoms to a 0.48mm diameter spot. Relay lenses 
for the hollow beams are not shown. Faraday pump and probe 
beams propagate along the x-axis to a balanced polarimeter. 
WP: Wollast on prism. Osc: Oscilloscope. Image of beam and 
cross sectional profile are shown. 
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Figure [T] shows the schematic layout of the experiment. 
The hollow beam is relayed to intersect itself by an 8/ 
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imaging relay, as described in Ref. [19]. Helmholtz coils 
on all three axes control the magnetic field. The hollow 
beams for our trap are formed by modifying the wave- 
front phase of a Gaussian beam with a reflective spa- 
tial light modulator (SLM). SLMs have found increas- 
ing value in cold atom manipulation experiments be- 
cause of their ability to control trap parameters in a 
programmable manner and to produce traps with non- 
trivial intensity profiles [H [201 EH El [^. The ap- 
plied phase for the hollow beams used here has a profile 
^(p, (j)) = ncj) -\- /A/(7rp^), where p and ^ are cylindrical 
coordinates and n is an integer. The second term is a lens 
function of focal length / ^ 200mm to focus the beam 
of wavelength A onto the atom sample. For high charge 
number beams (n > 4), we usually operate the trap a few 
centimeters away from the focal plane, where aberrations 
are reduced and the peak intensity is maximum [18]. 

The light for the hollow beam is derived from a tunable 
extended cavity diode laser. It is amplified to 400 mW 
by a tapered amplifier, 200 mW of which is coupled into 
polarization-maintaining fiber. Residual resonant light 
from amplified spontaneous emission is filtered out by a 
heated vapor cell. The fiber output is collimated to a 
1/e^ waist of 1.71 mm, and modified by the SLM (Boul- 
der Nonlinear Systems), which has ~ 90% diffraction ef- 
ficiency. The SLM has been calibrated at the pixel level 
to correct for wavefront distortion intrinsic to the SLM. 
An image of the beam is shown in the inset to Fig. [l] 
Our choice of n = 8 is driven by the practical consider- 
ations of field- free confinement and large trap size, but 
these trap parameters can be adjusted with the SLM. 

Our experiment begins with cold ^^Rb atoms derived 
from a magneto-optical trap (MOT). We confine ^ 10^ 
atoms in a ~ 500/im diameter (1/e^) cloud. The atoms 
are further cooled in a 10 ms long molasses stage to 
^ 10/iK, after which all MOT-related beams are extin- 
guished. The hollow beam trap is on throughout the 
MOT loading, but can be switched off by an acoustoop- 
tic modulator. The Faraday spectroscopy is performed 
by similar technique as in Ref. [2 . To perform these 
measurements, a pair of laser beams is used along the 
X-axis (Fig.jl]). The atoms are optically pumped into the 
F = 3, mi? = 3 stretched state by a 20 /is cr+ pulse 
connecting F = 3 ^ = 3. This beam has 1/e^ 
waist of 6.0 mm and has a peak intensity of ~ 31 sat ^ 
where I sat is 1.6 mW/cm^. This beam is retroreflected 
to prevent unidirectional momentum kicks, and a small 
amount of repumper light 0.04/sat) is added during 
this pulse to keep the atoms in the F = 3 hyperfine 
ground state. When this light is extinguished, the atoms 
begin precessing freely at the Larmor precession fre- 
quency ujl = grPsB/h^ where is the gyromagnetic ra- 
tio, /iB is the Bohr magneton, and B is the magnetic field. 
For ^^Rb, grps/h = 466.7415 kHz/Gauss [24]. A lin- 
early polarized probe beam at a detuning Ap = 27r x 2.5 
GHz with c^20mW and 1/e^ waist Up = 6.0mm passes 
through the atom cloud to a simple polarimeter consist- 
ing of a Wollaston prism that splits the probe beam into 



two orthogonal polarization states that are detected by a 
balanced photodetector. For these parameters, the pho- 
ton scattering time from the probe beam is calculated to 
be Tp ^2ms. The MOT region is imaged onto a pinhole 
along the axis of the probe beam so that only the por- 
tion of the probe that interacts with the confined atoms 
reaches the detector. 

The hollow beam trap prevents the atoms from falling 
away from the interaction region during the probing pro- 
cess. The beam has 150 mW total power at the trap. 
We use a detuning ^ 25GHz (=0.05 nm) above the 
F = 3 ^ = 4 transition. At the MOT, the hol- 
low beam has a diameter of 0.48 mm, measured between 
maxima, and the peak intensity is 8.2 x 10^ mW/cm^ 
for a trap depth U ^ 2hT ^ 3000F^, where Er is the 
recoil energy for ^^Rb. The gravitational potential en- 
ergy across this trap is HT/G. For these parameters, 
the peak scattering rate from the trapping beams would 
be 7^ = 1/rt ^ 27r x 3kHz, but is reduced from this 
value by being trapped in the dark. Although we do not 
measure this value, we establish an upper bound to be 
7t<27r X 200Hz. 
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FIG. 2: (Color online) Faraday signals for a) trapped sample; 
and b) untrapped sample. The untrapped atoms fall away 
from the probe region within 20 ms, while the trapped atoms 
remain with 150 ms time constant. Insets to (a) are expanded 
views of the raw data. 

Each optical pumping event initiates the Larmor pre- 
cession. Figure|2^ shows 64 averages of 200 optical pump- 
ing cycles spaced 2 ms apart in the presence of the hollow 
beam trap and a bias magnetic field of ~ 100 mG along 
the z-axis. A single Larmor precession signal is shown 
in the lower inset to Fig. [2^. The envelope over all Lar- 
mor precession signals decays with a 1/e time constant 
of ~ 150 ms. This decay is due primarily to the steady 
heating that occurs during each optical pumping cycle. 
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which graduahy boils atoms out of the trap. In contrast, 
Fig. (2)3 shows the signals without the hollow beam trap 
present. In this case, the atoms fall completely out of the 
probe beam detection window within 25 ms, with a 1/e 
decay time of 13 ms, similar to that reported in Ref. [2]. 
The signals in Fig. [2] are recorded immediately following 
the molasses phase of the MOT loading cycle. Over the 
first few pumping cycles, the envelope of the individual 
precession signals in Fig. (2)3 changes dramatically due to 
residual eddy currents in the vacuum chamber. Holding 
the atoms in an optical trap allows measurements to be 
performed after eddy currents have subsided, while also 
substantially increasing both the measurement window 
and the overall duty cycle. 

For our parameters, each independent Larmor preces- 
sion signal dephases with a submillisecond 1/e decay 
time. This dephasing occurs from several factors, in- 
cluding spatial gradients and photon scattering from the 
trap and probe beams. Nonlinear Hamiltonian terms can 
also shorten the decay time of the signal, as described in 
Ref. [12 . These nonlinear terms depend on the angle 
between the polarization of the probe laser and the mag- 
netic field. When the relative angle is ~ 54°, the effects 
of these terms are eliminated. For this work, we operated 
at this relative orientation so that the dephasing occurs 
primarily through photon scattering. From Fig. |2j we 
find that the untrapped signals decay with a 1/e time of 
?^0.7ms. For the samples trapped in the hollow beam, 
we observe a slight reduction in the decay time to ^0.5 
ms. Thus we have an upper bound for 7t < 27r x 200Hz. 
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FIG. 3: (Color online) a) Polarimeter output for single- 
shot data (top) and averaged data (bottom), b) Fast-Fourier 
transform (FFT) of data in (a). Fits to a Lorentzian profile 
are shown as solid lines. Single-shot (circles) and averaged 
(squares) data are shown. 

In a gradient-free, static magnetic field, the voltage 
output of the polarimeter is an exponentially-damped si- 
nusoid, V{t) = Aexp(-t/r)sin(27ri^i,t + 0), where A is the 
initial amplitude, r is the 1/e decay time, col = '^tt x 



is the Larmor frequency, and is a phase. To deter- 
mine z^L, the averaged data in each 2 ms probing window 
(Fig. [3^) are Fourier transformed (Fig.|3]3). We fit these 
transforms to a Lorentzian, the center of which is ul. 
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FIG. 4: (Color online) a) Larmor precession frequency as a 
function of time with various levels of compensation. Dashed: 
no compensation; dotted: compensation of eddy currents, and 
solid: full compensation. Inset to (a) shows magnified view 
of compensation. Fluctuations in our Larmor precession are 
dominated by uncompensated harmonics of the AC line, b) 
Magnetic field spectrum with no compensation (dotted), and 
60Hz compensation (solid). 

In Fig. [4^, we plot VL{t) over 200 pumping cycles (64 
averages) spaced 1 ms apart (dashed line). The sig- 
nal displays two dominant sources of time-dependence. 
First, the exponential decay occurs from the metal vac- 
uum chamber, which develops eddy currents when the 
MOT coils are extinguished. Due to the symmetry of 
our chamber, the eddy currents are along the axis of 
the MOT coils (z). The bias field of ~ 50kHz for 
these measurements is also on the 2:-axis so that the eddy 
current field adds linearly to the bias field. The second 
source of time-dependent behavior is ambient AC mag- 
netic fields in the room arising from power supply trans- 
formers, power strips, etc. We note that our experiment 
is triggered off the AC power line. We found that this 
field is also primarily along the z-axis, because the ampli- 
tude of the oscillation signal is independent of this bias 
field. An orthogonal component would add in quadrature 
and cause the amplitude to vary with the bias field. Ad- 
ditionally, an orthogonal, oscillating magnetic field com- 
ponent added in quadrature would show up at twice its 
oscillation frequency. Since the Larmor frequencies re- 
trieved from Faraday spectroscopy determine the scalar 
magnetic field, full vector information is not acquired in 
a single shot, but can be acquired through multiple mea- 
surements [27 . Some information about magnetic field 
orientation can be obtained directly from the polarimeter 
signal (e.g. there is no spin precession if the B field is 
parallel to the optical pumping axis), but that effect is 
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outside the scope of this work. 

For many apphcations, control over the magnetic field 
is required to sub-mG levels, especially those involv- 
ing Raman transitions between magnetically sensitive 
states [25l [26l [27l [28j. As a simple application of the 
long measurement time capability, we demonstrate com- 
pensation of these time varying fields. We first compen- 
sate the effects of eddy currents, which produce an ex- 
ponentially decaying magnetic field at the atom sample. 
This field decays with a 1/e time of ~20 ms (Fig. |4|l). 
For a given MOT coil current setting, the eddy current 
amplitude is constant. We produce an opposing time- 
varying field fiux by using a voltage-controlled current 
source (Kepco ATE15-15M). This current passes through 
a 20-turn Helmholtz pair of diameter 20 cm, width 2.5 
cm, and separation 11.4 cm oriented along the MOT coil 
axis. The appropriate time variation is done by low- 
pass filtering of a step function whose amplitude is ad- 
justed for optimum compensation. The result is shown 
in Fig. ^ (dotted line). Although this source of time 
variation is not canceled perfectly, the field beyond 25 
ms is constant to within the 60 Hz field amplitude. 

The ambient AC magnetic fields are primarily due to 
60 Hz power line sources. By triggering our experiment 
from the power line, this source of magnetic field varia- 
tion is reproducible and can be compensated. Without 
this triggering, the variations of a few mG observed in 
Fig. [4^ would lead to significant shot-to-shot fiuctuations 
of the field measurements. We produce an opposing field 
by adding a 60 Hz sinusoidally varying current to the bias 
coils. The current amplitude and phase are adjusted for 
optimum compensation. The result with all compensa- 
tions applied is shown in the inset to Fig. [4^. The signal 
remains constant to within a standard deviation of 110 
Hz (230 /iG). Most of this residual field is due to higher 
AC line harmonics; in Fig. |4]3, we show frequency spec- 
trum of the magnetic field, which clearly shows higher 
harmonics at 180, 300, and 420 Hz. We suppressed the 
60 Hz component by a factor of 20. With appropriate 
signal processing, the field measurements in our setup 
could be made in real time (with single-shot measure- 
ments as in Fig. [3| and be used as feedback control with 
a bandwidth determined by the Helmholtz compensation 
coils. 

Another way to visualize the time dependent signals, 
shown in Fig.|5^, is by converting the ID data set of Fig.|2] 
to a 2D matrix. Each successive column contains the 
Larmor precession signal for subsequent triggers. This 
exposes time variations in an easily identifiable way with 
no FFT analysis. We show these images for the magnetic 
fields with no compensation, 60 Hz compensation and 
full compensation. A constant magnetic field shows up 
as a series of horizontal lines whose spacing is inversely 
proportional to ljl (Fig. [sJ)). 

Because of the uncompensated field variations in our 
lab, our measurement uncertainty is dominated by sys- 
tematic errors. To differentiate the systematic error from 
the random error, we measure the Larmor precession fre- 
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FIG. 5: 2D image of data that exposes qualitative magnetic 
field features without FFT processing, (a) uncompensated, 
(b) 60 Hz compensation, c) full compensation. 
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FIG. 6: (Color online) Comparison of 64 data runs (dashed) 
to a simulated signal with white noise added to have similar 
SNR as the experiment (black). For clarity, the signals are 
offset from each other. 



quency in a 1 ms measurement window for 64 indepen- 
dent loading cycles at a trapping time of T = 100 ms. 
These scans are recorded at a 1 Hz rate. The result is 
shown in Fig. [6j There is a long term drift in our lab 
on the order of several seconds. Our time-domain signals 
have a signal-to-noise ratio (SNR) of ^15 (Fig. [3|i). For 
comparison, we simulated the expected Larmor signals 
for exponentially-damped sinusoids with the same SNR 
that had additive white noise (Fig. [6|. For the experi- 
mentally measured case, we found a standard deviation, 
or single-shot error, in each 1 ms optical pumping cy- 
cle of SuL ^45Hz, or SB ^ 100/iG (=10 nT), and for the 
simulated case, the error was Sul ^16 Hz, or SB ^ 30/iG 
(=3 nT). This discrepancy is likely due to other sources, 
such as unwanted variations in the MOT coil current. 

Our hollow beam traps are initially loaded with 
N ^ 10^ atoms. The shot-noise-limited magnetic 
field measurement error due to atom number is SB 



(h/gjiB) (l/V^^Vrl^), where r is the spin-coherence time 
and Tm is the measurement time [1 . Because we are 
measuring a rapidly varying field, Tm = 1— 2ms, limit- 
ing SB ^ 2/iG (=200 pT) in each optical pumping cycle. 
After T = 400 ms of trapping time, when there are only 
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^ 10^ atoms remaining, this increases to ^ 6/iG (=600 
pT). Our measured values are above the shot noise hmit 
due to the simple photodetection circuit we used and to 
incomplete optical pumping, which effectively reduces N. 

For static magnetic fields, each measurement cycle 
through the total trap time can be averaged, effectively 
increasing to several hundred milliseconds and greatly 
increasing the shot-noise-limited sensitivity. Likewise, r 
can be increased by using larger detunings for the probe 
and trapping beams. These blue-detuned traps are capa- 
ble of capturing large enough atom numbers that mea- 
surements in the low pT range or better should be pos- 
sible in a single MOT loading cycle over the entire mea- 
surement window. 
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FIG. 7: (Color online) Simulations of atom number remain- 
ing for different degrees of optical pumping. With no optical 
pumping or probe beams (longest- lived curve) the atoms still 
scattered ^1 photon/ms from the trap beams. A curve fit 
is shown for the case of scattering 7 photons/ms. This time 
constant of 160 ms agrees well with our experimental value of 
150 ms (Fig.[^. 

For any trap depth, there is a trade-off between SNR 
and the number of possible field measurements allowed 
before the signal decays. SNR improves by increasing 
the number of atoms that are optically pumped or by 
decreasing the pump detuning [5 , but these approaches 
also boil the atoms out of the trap more quickly. For most 
of our results presented here, we only weakly pumped the 
atoms to reduce heating and to increase the number of 
optical pumping cycles we could achieve. In general, the 
dominant heating will occur from the 20/is optical pump- 
ing phase of each cycle, during which several photons are 
scattered. As a rough estimate, the timescale for signal 
decay should be on the order of the time required for the 
average atom energy to equal the trap depth. This will 
occur after a time Tijou = U/ {jtotEr), where U is the trap 
depth, 7tot is the total scattering rate (including probe, 
trap, and optical pumping beams), and Er is the recoil 
energy. For our trap ofU ^ 3000£^r, and assuming ^ 10 
scattered photons every 2ms optical pumping cycle, this 
gives ~ 300 pumping cycles before the atoms are boiled 
away. 

To examine this boiling process more accurately, we 
perform Monte Carlo simulations of the atom dynam- 
ics within our trap for different total scattering rates. 
Within each time step, the atom's momentum is changed 



with a probability determined by the local scattering 
rate for the probe and trap beams, as calculated by the 
Kramers- Heisenberg formula [29]. We performed these 
simulations for different optical pumping rates. In Fig. [T] 
we plot the number of atoms remaining as a function of 
time for varying degrees of scattering rates. For the case 
of 7 photons scattered every millisecond, we find an expo- 
nential decay of ^160 ms, which is close to our observed 
value of Ti)oii = 150ms. By turning off the probe and 
optical pumping beams in the simulation, we find that 
the trap beam scattering rate is 7^ ~ 27r x lOOHz which 
agrees with our upper bound of 7^ < 27r x 200Hz from 
the Faraday decay time. 

Within our measurement error, we observed no effect 
of the trapping light on the Larmor frequency. Optically- 
induced Zeeman shifts that occur with elliptically polar- 
ized light [71 Ull [30] should be small, because the trap 
beam polarizations are linear and because of the low field 
confinement. Furthermore, any vector light shifts from 
the trap beams, confined to the x — y plane, would add 
in quadrature to our applied magnetic field along re- 
ducing the effect on col [H^- We are currently studying 
the effects of trap geometry on the Larmor precession 
signals. 
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FIG. 8: (Color online) Larmor precession signals with and 
without the trap light at two different relative orientations 
of the laser polarization with respect to the magnetic field. 
These were done at a) degrees (revivals maximized) and b) 
54 degrees (revivals suppressed). 

The dominant source of nonlinear effects due to the 
laser fields is the probe light. As discussed in Ref. [ 12] . 
the tensor component to the light shift adds a nonlin- 
ear term to the spin Hamiltonian, whose magnitude is 
dependent on the angle between the laser polarization 
and the magnetic field. This Hamiltonian plays an im- 
portant role in studies of quantum chaos and is use- 
ful for both nondestructive quantum state preparation 
and measurement [12]. For sufficiently large magnetic 
fields, the nonlinearity vanishes when the relative angle 
is ^=arctan(A/2) ^ 54°, but is maximized for ^ = 0. We 
have verified that these nonlinear spin dynamics, which 
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manifest themselves as revivals of the Faraday oscilla- 
tion signal, can still be observed in these hollow beam 
traps. In Fig. |8j we show the Larmor precession signals 
for ^ = and = 54° both with the trap on continu- 
ously and with the trap switched off immediately prior 
to the optical pumping pulse. Thus the high duty cycle 
of the technique presented here may be of use for rapidly 
testing quantum state preparation procedures employing 
this nonlinearity. 

We have demonstrated Faraday spectroscopy with high 
repetition rate, long measurement time, and submillime- 
ter spatial resolution in a dark hollow beam optical trap. 



We used high- charge- number hollow laser beams to pro- 
vide box-like confinement with near resonant light and 
low laser power. These traps can be sufficiently deep 
that several hundred Faraday measurements are possible 
before atoms are heated over the confining potential. We 
demonstrated a continuous magnetic field measurement 
over a period of 400 ms which enabled us to measure and 
compensate for time- varying magnetic fields. This work 
was funded by the Office of Naval Research and by the 
Defense Advanced Research Projects Agency. 
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